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00 ' The ionization of atomic electrons by scattering of neutrinos is revisited. This process is the one 

OO ■ studied in the experimental searches for a neutrino magnetic moment using germanium detectors. 

Current experiments are sensitive to the ionization energy comparable with the atomic energies, 
and the effects of the electron binding should be taken into account. We find that the so-called 
stepping approximation to the neutrino-impact ionization is in fact exact in the semiclassical limit 
^ . and also that the deviations from this approximation are very small already for the lowest bound 

H : 

, Coulomb states. We also consider the effects of electron-electron correlations and argue that the 

resulting corrections to the ionization of independent electrons are quite small. In particular we 
estimate that in germanium these are at a one percent level at the energy transfer down to a 
fraction of keV. Exact sum rules are also presented as well as analytical results for a few lowest 
hydrogen-like states. 
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I. INTRODUCTION 



The neutrino magnetic moments (NMM) expected in the Standard Model are very small 
and proportional to the neutrino masses /i^ ~ 3 x 10^^^ /i^ eV) with /i^ = e/2m 

being the electron Bohr magneton, and m is the electron mass. Thus any larger value of 
/ij, can arise only from physics beyond the Standard Model (a recent review of this subject 
can be found in Ref. i2||). Current direct experimental searches 3-5| for a magnetic moment 
of the electron (anti)neutrinos from reactors have lowered the upper limit on fi,y down to 
fii, < 3.2 X 10~^^ [5]. These ultra low background experiments use germanium crystal 
detectors exposed to the neutrino flux from a reactor and search for scattering events by 
measuring the energy T deposited by the neutrino scattering in the detector. The sensitivity 
of such a search to NMM crucially depends on lowering the threshold for the energy transfer 
T, due to the enhancement of the magnetic scattering relative to the standard electroweak 
one at low T. Namely, the differential cross section da/dT is given by the incoherent sum 
of the magnetic and the standard cross section, and for the scattering on free electrons the 



NMM contribution is given by the formula 



la 



dT " \T Ey J m? \fJ,B J \T 

where Ey is the energy of the incident neutrino, and displays a 1/T enhancement at low 
energy transfer. The standard electroweak contribution is constant in T aX E^ ^ T: 



daEW GpTTl , , A ■ 2a I C • 4/) \ 



dT 271 



1 + 



T 
E'y 



5 X 10"^^ cmVkeV. (2) 



In what follows we refer to these two types of contribution to the scattering as, respectively, 
the magnetic and the weak. 

The current experiments have reached threshold values of T as low as few keV and are 
likely to further improve the sensitivity to low energy deposition in the detector. At low 
energies however one can expect a modification of the free-electron formulas ([1]) and ([2]) 
due to the binding of electrons in the germanium atoms, where e.g. the energy of the Ka 
line, 9.89 keV, indicates that at least some of the atomic binding energies are comparable to 
the already relevant to the experiment values of T. Thus a proper treatment of the atomic 
effects in neutrino scattering is necessary and important for the analysis of the current 
and even more of the future data with a still lower threshold. Furthermore, there is no 
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known means of independently calibrating experimentally the response of atomic systems, 
such as the germanium, to the scattering due to the interactions relevant for the neutrino 
experiments. Therefore one has to rely on a pure theoretical analysis in interpreting the 
neutrino data. This problem had been addressed numerically in the past ^, 9| and the 
interest to the problem was renewed in several recent papers, which however are ridden by 
a 'trial and error' approach. The early claim 10|] of a significant enhancement of the NMM 



contribution by the atomic effects has been later disproved 



11 



12| and it was argued 



that the modification of the formulas ([T]) and ([2]) by the atomic binding effects is insignificant 



down to very low values of T. It has been subsequently pointed out [13| that the analysis 
ot Rat. Q i. generally „„.l.dated ^n,ulti-elect.o„ =y=tem, Muding atoms „.th Z > 1. 
Furthermore, the analysis of Ref. [llj is also generally invalidated by singularities of the 
relevant correlation function in the complex plane of momentum transfer, so that the claimed 
behavior of the cross section at low T applies only in the semiclassical limit, although, as 
will be shown here, it gives a very good approximation to the actual behavior for an electron 
bound by a Coulomb potential. 

In this paper we revisit the subject of neutrino scattering on atoms at low energy transfer. 
We aim at describing this process at T in the range of few keV and lower, so that the motion 
of the electrons is considered as strictly nonrelativistic. Also in this range the energy of the 
dominant part of the incident neutrinos from the reactor is much larger than T and we thus 
neglect any terms whose relative value is proportional to T / (in particular, in this range 
one can neglect the 1/Ey term in Eq.([l]) in comparison with 1/T). Furthermore any recoil 
of the germanium atom as a whole results in an energy transfer less than 2El/MGei which at 
the typical reactor neutrino energy is well below the considered here keV range of the energy 
transfer. Thus we formally set the mass of the atomic nucleus to infinity and neglect any 
recoil by the atom as a whole. In particular, under these conditions the interaction of the 
neutrino with the nucleus can be entirely neglected, and only the scattering on the atomic 
electrons is to be considered. 

We find that in the scattering on realistic atoms, such as germanium, the so-called step- 
ping approximation works with a very good accuracy. The stepping approach, introduced 
in Ref. [8i] from an interpretation of numerical data, treats the process as scattering on 
individual independent electrons occupying atomic orbitals and suggests that the cross sec- 
tion follows the free-electron behavior in Eqs.([T]) and ([2]) down to T equal to the ionization 
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threshold for the orbital, and that below that energy the electron on the corresponding or- 
bital is 'deactivated' thus producing a sharp 'step' in the dependence of the cross section on 
T. In the present paper, we consider general relations for the discussed scattering on atomic 
systems in Sec.[ITland present in Appendix |X] sum rules for the theoretical objects involved 
in the calculations 18| . In Sec. IIIII we prove that for the scattering on individual electrons 



the stepping approximation becomes exact in the semiclassical limit, so that its applicability 
is improved with the principal number n of the atomic orbital. We also find by an explicit 
calculation (Appendix [B]) for a hydrogen-like ground state, i.e. at n = 1, that the deviation 
from the stepping behavior is less than 5% at the worst point, where the energy transfer T is 
exactly at the threshold. The accuracy of the approach based on considering the scattering 
on individual electrons is limited by the existence of the electron-electron correlations in 
the process. We consider the correction introduced by these correlations in Sec. |IV] and, 
in Sec. |Vl apply the derived formula to an estimate of the effect for germanium, using the 
Thomas-Fermi model. We find that the correlation correction grows at smaller T but is still 
small, of order of a few percent, for T in the range of a few hundred eV. We thus argue that 
the stepping approach describes the scattering cross section with a sufficient for practical 
purposes accuracy, and that it can be applied to the analysis of the present and future data 
of searches for NMM with germanium detectors down to the values of the energy deposition 
T ~ 0.3 keV. 



II. GENERAL FORMULAS FOR NEUTRINO SCATTERING ON ATOMIC 
ELECTRONS 

In this section we briefiy recapitulate the general expressions and introduce the relevant 
atomic objects for the neutrino scattering on atomic electrons. We start with the magnetic 
process and then also apply a similar treatment to the standard weak part of the cross 
section. 

The kinematics of the scattering of a neutrino on atomic electrons is generally char- 
acterized by the components of the four-momentum transfer, the energy transfer T and 
the spatial momentum transfer g, from the neutrino to the electrons with two rotationally 
invariant variables being T and = q^. At small T the electrons can be treated nonrela- 
tivistically both in the initial and the final state, so that the process is that of scattering of 
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an NMM in the electromagnetic field A = {Aq,A) of the electrons: Ao^q) = y/iira piq)/q^, 



A(q) = y/iTia j (q) / , where p{q) and j{q) are the Fourier transforms of the electron number 
density and current density operators, respectively, 

z 



p{q) = ^exp{iq-fa) , 



a=l 



a=l 



exp(zg ■ rajTT^ + 77^ exp(zg ■ Ta) 



(3) 



(4) 



and the sums run over the positions of all the Z electrons in the atom. 

In this limit the expression for the double differential cross section is given by 



13| 



dT dq^ 



(-?) 



^(T,g2)+(l-^)/2(T,g^) 



(5) 



where S{T,q ), also known as the dynamical structure factor [15|, and R{T,q ) are 

SiT, q') = Y^5{T-E^ + Eo) |(n|p(g)|0)|^ , 

n 

i?(T, q') = J25{T-En + Eo) \{n\umO)f , 



(6) 



(7) 



with jj_ being the j component perpendicular to q and parallel to the scattering plane, which 
is formed by the incident and final neutrino momenta. The sums in Eqs. ([6]) and ([7]) run 
over all the states |n) with energies En of the electron system, with |0) being the initial 
state. 

Clearly, the factors S(T, q^) and -R(T, q^) are related to respectively the density- density 
and current-current Green's functions 

1 



L{T,q') = J2 



\{n\p{mr 

T-En + Eo-te 

\{n\j±{m\" 
T-En + Eo-ie 







T-H + Eo-ie 
1 

T -H + Eo-ie 







as 



S{T,q^) = - ImF(T, q^ 

TT 



R{T,q' 



TT 



ImL(T, g2) 



(9) 

(10) 
(11) 
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with H being the Hamiltonian for the system of electrons. For small values of g, in particular, 
such that g ~ T, only the lowest-order non-zero terms of the expansion of Eqs. ffTOj) and f lTT ) 
in powers of are of relevance (the so-called dipole approximation). In this case, one has [13 1 



R{T,q') = ^S{T,q'). (12) 

Taking into account Eq. (|T2l) . the experimentally measured singe-differential inclusive 
cross section is, to a good approximation, given by (see e.g. in Refs. [ll|,Q) 

^ = 4vra/.,y^^ 5(T,g)— . (13) 

The standard electroweak contribution to the cross section can be similarly expressed in 
terms of the same factor S{T, q^) [lll| as 

^ = ^ (1 + 4 sin^ + 8 sin^ 9^) S{T, q') dq' , , (14) 

where the factor S{T, q^) is integrated over with a unit weight, rather than as in 
Eq.dlS]). 

The kinematical limits for q^ in an actual neutrino scattering are explicitly indicated in 
Eqs.( IT3|) and (fT4|) . At large Ey, typical for the reactor neutrinos, the upper limit can in fact 
be extended to infinity, since in the discussed here nonrelativistic limit the range of momenta 
~ is indistinguishable from infinity. The lower limit can be shifted to = 0, since the 
contribution of the region of q^ < can be expressed in terms of the photoelectric cross 
section jll| and is negligibly small (at the level of below one percent in the considered range 
of T). For this reason we henceforth discuss the momentum-transfer integrals in Eqs. ( !T3|) 
and ( 1T4|) running from g^ = to g^ = oo: 

h{T) = / 5(r, g2) 4 , and h{T) = / S{T, q') dq' . (15) 

Jo q Jo 

For a free electron, which is initially at rest, the density- density correlator is the free 
particle Green's function 

2m 

so that the dynamical structure factor is given by S(^pE)(T,q'^) = 5{T — g^/2m), and the 
discussed here integrals are in the free-electron limit as follows: 

= [ S^,j,){T, q')^ = ^, = [ S^.E^iT, q^) dq^ = 2m. (17) 



F(^^)(T,g^)=(r-|--.e) (16) 



Clearly, these expressions, when used in the formulas ( !T3|) and ( IT^ . result in the free-electron 
cross section in Eqs. ([1]) and (EJ. 

III. SCATTERING ON ONE BOUND ELECTRON 

The binding effects generally deform the density- density Green's function, so that both 
the integrals f lTSj) are somewhat modified. Namely, the binding effects spread the free- 
electron (5-peak in the dynamical structure function at = 2mT and also shift it by the 
scale of characteristic electron momenta in the bound state. However it turns out that 
the free electron expressions are quite robust in the sense that in realistic systems the 
modification of the integrals relative to their free-electron limit, are quite small. As a formal 
statement, we will show in the Appendix |X] that when the function F(T, g^) is analytically 
continued in complex plane of the free-electron expressions are valid for the integrals over 

extending from — oo to +00, and in the case of the integral, similar to I2 i.e. with the 
weight this property also holds for scattering on multi-electron atomic systems, while 
for that with the weight it generally holds only for the scattering on one electron, or on 
independent electrons. Clearly the latter integrals over the full axis of differ from those of 
physical interest in Eq. (fT5|) by the contribution of negative g^, which although numerically 
small even at low T, still makes the scattering on bound electrons different from that on 
free electrons. 

In this section we consider the scattering on just one electron. The Hamiltonian for the 
electron has the form H = /2m + V(r), and the density- density Green's function from 
Eq-dHD can be written as 



where the infinitesimal shift T — t- T — ze is implied. 

Clearly, a nontrivial behavior of the latter expression in Eq.( |T8l) is generated by the 
presence of the operator {p ■ q) in the denominator, and the fact that it does not commute 
with the Hamiltonian H. Thus an analytical calculation of the Green's function as well 
as the dynamical structure factor is feasible in only few specific problems. In Appendix [B] 



F(T,g2) 



(0 I e'''^-'' [T - H{p, f) + Eo]-^ e''^-' \ O) 
{0\[T-H{p + q,r) + Eo]-'\0) 




(18) 
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we present such a calculation for ionization from the Is, 2s and 2p hydrogen-like states. 
In particular, we find that the deviation of the discussed integrals f|T5|) from their free 
values are very small: the largest deviation is exactly at the ionization threshold, where, for 
instance, each of the Is integrals is equal to the free-electron value multiplied by the factor 



;i -7e-V3) ^0.957 [l9|. 



2m m J \ 2m m 



(19) 



The problem of calculating the integrals ( ITSjl however can be solved in the semiclassical 
limit, where one can neglect the noncommutativity of the momentum p with the Hamil- 
tonian, and rather treat this operator as a number vector. Taking also into account that 
{H — Eq) |0) = 0, one can then readily average the latter expression in Eq. ([TK|) over the 
directions of q and find the formula for the dynamical structure factor: 

2pq 

where p = \p\ and 6 is the standard Heaviside step function. The expression in Eq. (fT9|) 
is nonzero only in the range of q satisfying the condition —pq/m < T — q^ /2m < pq/m, 
i.e. between the (positive) roots of the binomials in the arguments of the step functions: 
qmin = a/ 2m T + p"^ — p and qmax = \/2mT + p^ + p. One can notice that the previously 
mentioned 'spread and shift' of the peak in the dynamical structure function in this limit 
corresponds to a flat pedestal between g^m and qmax- The calculation of the integrals ( ITSjl 
with the expression f lT9|) is straightforward, and yields the free-electron expressions ([T7|l for 
the discussed here integrals in the semiclassical (WKB) limit j^oj]: 

Jr''^) = i, Jf^^) = 2m. (20) 

The difference from the pure free-electron case however is in the range of the energy trans- 
fer T. Namely, the expressions (l20ll are applicable in this case only above the ionization 
threshold, i.e. at T > |£'o|- Below the threshold the electron becomes 'inactive'. 

It is instructive to point out that the validity of the result in Eq. fl20|) is based on the 
semiclassical approximation and is not directly related to the value of the energy T. In 
particular, for a Coulomb interaction the WKB approximation is applicable at energy near 
the threshold [l^. For T exactly at the threshold, T = — -Eq? the criterion for applicability 
of the semiclassical approach in terms of the force F = \F\ = \ [p,H] \ acting on the electron 



and the momentum p of the electron is that iJ] the ratio of the characteristic values mF/p 



.3 



is small. For the excitation of a state with the principal number n this ratio behaves para- 



metrically as 1/n 2l| . Thus the applicability of a semiclassical treatment of the ionization 
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near the threshold improves for initial states with large n. As previously mentioned, the 
modification of the integrals f lTSj) by the binding is already less than 5% for = 1, so that 
we fully expect this deviation to be smaller for the higher states, and even smaller at larger 
values of T above the threshold due to the approach to the free-electron behavior at T ^ Eq. 

We believe that the latter conclusion explains the so-called stepping behavior observed 
empirically jsj] in the results of numerical calculations. Namely the calculated cross section 
da/dT for ionization of an electron from an atomic orbital follows the free-electron depen- 
dence on T all the way down to the threshold for the corresponding orbital with a very small, 
at most a few percent, deviation. This observation lead the authors of Ref. to suggest 
the stepping approximation for the ratio of the atomic cross section (per target electron) to 
the free-electron one: 

fm^J^^j^^-^En-nT-m), (21) 

where the sum runs over the atomic orbitals with the binding energies Ei and the filling 
numbers n^. Clearly, the factor /(T) simply counts the fraction of 'active' electrons at the 
energy T, i.e. those for which the ionization is kinematically possible. For this reason we 
refer to f{t) as an activation factor. We conclude here that the stepping approximation is 
indeed justified with a high accuracy in the approximation of the scattering on independent 
electrons, i.e. if one neglects the two-electron correlations induced by the interference of 
terms in the operator p{q) in Eq.Q corresponding to different electrons. In the next section 
we estimate the effect of such an interference and find that the resulting corrections are 
small, at least in atoms with large Z, such as the germanium. 



IV. TWO-ELECTRON CORRELATION 



In this section we discuss the correction arising from a correlation between two electrons. 
We consider the energy T and hence the relevant momentum transfer q as large in comparison 
with the atomic scale. In this way we estimate the relevant parameter for the significance 
of the correlation effect. 

We start with consider an isolated system of two electrons interacting among themselves 
through the Coulomb potential V{r). The Hamiltonian for this system thus has the form 

p2 2 

H=— + ^ + V{r) , (22) 



where P = —id/dR and p = —id/df are, as usual, the momenta conjugate to respectively 
the center of mass coordinate R and the relative coordinate f. 

The spatial part of the wave function of the system factorizes into the product 0(-R) tp{r) 
with (f){R) while the spin part will be considered later. We consider here the system at rest, 
i.e. 0(-R) = const, since a boost to a momentum P does not change the cross section. 

The density- density Green's function then takes the form 

nr, q') = E ^ ^ ^ — ' (23) 







where the states |0), \n) and the energies Eq, En refer to the relative motion in the system 
with |0) standing for the initial state. Clearly, it is implied in Eg. (123 p that the corresponding 
matrix elements for the (trivial) dynamics of the system as a whole are already taken, which 
results in replacing in the energy denominator the excitation energy T by its value corrected 
for the recoil of the system as a whole: T — t- T — /Am. 

The cross terms between exp(zg* ■ r/2) and exp(— ■ r/2) in the expression fl2^ result in 
the previously discussed one-particle Green's function 



Fi(T,,=) = 2 (0 



2m m 



, (24) 



where the overall factor of 2 arises from the two identical (after averaging over the direction 
of q) cross terms, and physically is corresponding to the presence of two particles in the 
system. 

The discussed here contribution of the two-electron correlation arises from the diagonal 
terms, whose contribution is given by 

(0 |e*9^/2| ^\ |gigr/2| Q\ 

^ T-^-E^ + Eo 

where, again, the two terms give the same contribution after the averaging over the direction 
of q, which is accounted for by the factor of 2 in the latter expression. The factor rj is the 
symmetry factor for the spin part of the two-electron system: rj = —1 for the spin-singlet 
state of the pair and = +1 for the spin-triplet. The appearance of this factor can be 
explained as follows. The discussed correlation arises from the situation where an excitation 
of one electron by the operator p produces the same spatial wave function as an excitation 
of another one. In order for the wave functions to be identical the spin variables of the 
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two electrons should also be switched, which operation results in the factor rj. Clearly, 
no such factor arises in the one-particle term since the spin of both electrons simply 
'goes through'. It can be also mentioned that, naturally, the symmetry of the spatial wave 
function ^/'(r) is opposite to r]. 

One can notice that unlike in the one-particle contribution (Eq.( l24p ) where the momen- 
tum q flows in and out of the system, the correlation contribution in Eq. (l25|l corresponds 
to the net momentum q flowing into the system. Clearly, for non-interacting particles such 
contribution would vanish and the whole correlation effect arises only due to the interaction 
between the electrons, which interaction absorbs the momentum transfer. The term f l25p 
can be graphically represented as shown in Fig.l, where the system lines correspond to the 
propagation of the system in the potential V with the outer legs corresponding to the wave 
functions of the initial state in the momentum space {ip{p') \ and \ip{p)) and the line between 
the action of the operators exp{iq ■ r/2) corresponding to the Green's function. One can 
write in terms of these objects the expression for Fc{T, q^) as 

^^^^'^') = '/ 1^ / ^.r^f^')G{E;p'-y+l)m , (26) 

where G{E^ k', k) is the Green's function {E — H)~^ in the momentum representation at the 
energy E = T + Eq — q^ /Am. 
We shall consider separately the effect of the interaction in the wave functions and in the 




FIG. 1: Graphical representation of the two-electron correlation. The external legs correspond to 
the momentum-space wave function and the propagator is the Green's function in the potential V . 

Green's function. For the zeroth order Green's function 

P ^ {2^)H^^\k'-k) 
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(27) 



and the exact wave functions one finds 

i^Ty) = 2, f-pL — riPpim — . (28) 

^ ^ J (2'Ky T - - ^ - ^ + Ea 

Let us consider now g as a large parameter in comparison with the characteristic momenta 
Po in V'(p)) beyond which the wave function falls off. At such values of q the product 
ip*{p + (l)'ip{'p) carries a suppression in only one of the factors in two regions of p: one where 
p ^ Po and the other where \p+q\ ~ Po- Clearly, by shifting the integration variable p+q — )■ p 
one can readily see that both integration regions the contribution of the latter region is the 
same as of the first one, so that one evaluates the integral in Eq. fl28p by considering only 
the contribution of the region p q and taking it with a factor of two. Then the leading at 
large q expression for the function in Eq. fl28|) takes the form 



Fl'W).4„fS^ f |^,(p,.4„^Wl, (29) 

2m 2m 

The appearance of the wave function at the origin f = 0, 1^(0), in this expression implies that 
at large q the considered contribution to the correlation arises only for an S-wave relative 
motion within the electron pair, which thus have to be a spin-singlet, and therefore rj = —1. 

In fact for an for an 5- wave motion the momentum-space wave function ip{q) can be also 
expressed at large q in terms of the position-space wave function at the origin ipi^)- Indeed, 
in the S-wave the wave function is a function of r: ilj{r). At small r the Coulomb repulsion 
between the electrons dominates over all other interactions and the Schrodinger equation 
for '?/'(r) reads 

1 2 (T 

^" ^'(r) + -^{r) = E^{r). (30) 

m m r r 

By requiring the two singular at r — > as 1/r terms to match in this expression, one finds 
that the derivative of ■?/'(r) over r at the origin is expressed in terms of ip{0): 

^'(0) = ^^(0). (31) 

A finite derivative over r implies that the gradient Vtp{r) = ril)'{r)/r is singular at r = 
0, so that the asymptotic at large |g| behavior of the momentum space wave function is 
proportional to 1/|^^ with the coefficient determined by '?/''(0), which in turn is determined, 
according to Eq.l lHT]) . by ip(fi): 

47Tma tp{0) 

w 
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Using this relation in Eq.( l29|) one finds in the large limit 



16 TT a 



m|V^(0)p 



T — -3— 

2m 



(33) 



The latter expression is manifestly proportional to first power of the interaction between 
the electrons. Therefore a similar contribution can arise from the first order in the expansion 
of the Green's function in the interaction potential V. In this order one finds for the discussed 
correlation part of F{T, q^): 

<fp' f d^p ip*{p')V{—q + p' — p)'ilj{p) 



27] 



(27r)3 J (27r)3 



2m m 



/m m m ^ 



(34) 



where V{k) is the Fourier transform of the potential, so that for the Coulomb repulsion 
between the electrons 



V{k) 



Ana 

U2 ■ 



(35) 



Considering as before the limit of large q and thus neglecting p and p' in comparison with 
q, one readily finds that the result is again proportional to |?/'(0)p, so that the effect remains 
only in the 5*- wave (and hence 1] = —1): 

\m\' 



F^\T,q') 



I mT ,q'^'^p'^ ,p''^ 



— 57r a- 



(36) 



Collecting the formulas (155]) and (15^ together one finds the estimate of the two-electron 
correlation part of F{T, q^) in the limit of large T and q^: 



Fc{T,q' 



^7r«|^(0)p 



2 m 



qUT-^ 



2 

2m 



2m I 



(37) 



The corresponding two-electron correlation correction to the integrals for the neutrino scat- 
tering cross section, is then calculated by shifting T — )■ T — ie and considering only the 
contribution of the singularity at = 2m T: 



1 



ImF,(T,g2) 



dq' 



—An a 



\m\' 



^ j lmF,{T,q'')dq^ = 



(38) 



TV J q^ mT^ TV 

Notice, that the two-electron contribution to the integral relevant for the standard elec- 
troweak scattering vanishes in the discussed approximation due to a cancellation of the two 
terms in Eq.([3j 
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The discussed here calculation shows, as expected, that at large q the two-electron corre- 
lation arises only when the electrons are separated by a short distance. For this reason one 
can relax the assumption, we made in the beginning of this section, that the system of two 
electrons is in a free motion. Indeed the same result would apply in the situation, where 
the pair as a whole moves in a potential that is sufficiently smooth so that the 'tidal force' 
interaction with the rest of the atomic system does not overcome the Coulomb singularity 
of the repulsion between the electrons at distances of order 1/q. 



V. SCATTERING ON ATOMIC ELECTRONS IN GERMANIUM 

In considering the neutrino scattering on actual atoms one needs to evaluate the depen- 
dence of the number of active electrons on T and generally also evaluate the effect of the 
two-electron correlations. The energies of the inner K, L and M orbitals in the germanium 



atom are well known (see e.g. Ref. {q] and references therein) and provide the necessary 
data for a description of the neutrino scattering by the stepping formula fl21l) down to the 
values of the energy transfer T in the range of the binding of the M electrons, i.e. at 
T > \Em\ ~ 0.18 keV. The corresponding steps in the activation factor are shown in Fig. 2. 
It can be mentioned that if one applies formulas of the Appendix [B] to the onset of the K 
shell step, i.e. just above 10.9 keV, the difference from the shown in the plot step function 
would be practically invisible in the scale of Fig. 2. 

Our goal in this section is to estimate the effect of the two-electron correlations in the 
scattering on germanium. We shall estimate this effect by considering the atomic number 
Z as a large parameter and using the Thomas-Fermi model, which, in spite of its known 
shortcomings, appears to be appropriate for evaluating average bulk properties of atomic 
electrons at large Z, such as in the problem at hand. 

In the Thomas- Fermi model (see e.g. Ref. |lj]) the atomic electrons are described as a 
degenerate free electron gas in a master potential (j){r) filling the momentum space up to 
the zero Fermi energy, i.e. up to the momentum pQ{r) such that pl/2m — e0 = 0. The 
electron density n{r) = pI/IStt"^) then determines the potential 0(r) from the usual Gauss 
equation. In the discussed picture at an energy transfer T the ionization is possible only 
for the electrons whose energies in the potential are above — T, i.e. with momenta above 
j9T(r) with p'^/2m — ecj) = —T. The electrons with lower energy are inactive. Calculating 
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the density of the inactive electrons as pj./(37r^) and subtracting their total number from Z, 
one readily arrives at the formula for the activation factor, i.e. the effective fraction of the 
active electrons Zf.'g^jZ as a function of T: 



where is the Thomas-Fermi function, well known and tabulated, of the scaling variable 



and, finally, xo(T) is the point where the integrand becomes zero, i.e. corresponding to the 
radius beyond which all the electrons are active at the given energy T . The energy scale 
To in germanium (Z=32) evaluates to Tq ~ 3.1 keV. The Thomas-Fermi activation factor 
for germanium calculated from the formula (l39l) is shown by the dashed line in the plot of 
Fig. 2. One can see that in the shown energy range it reasonably approximates the stepping 
behavior of the atomic orbitals. The discussed statistical model is known to approximate 
the average bulk properties of the atomic electrons with a relative accuracy 0(Z~^/'^) and as 
long as the essential distances r satisfy the condition Z^^ ^ mar ^ 1, which condition in 
terms of the scaling variable x reads as Z'"^!^ <^ x <^ Z^l"^ . In terms of the formula (|39l) for 
the number of active electrons, the lower bound on the applicability of the model is formally 
broken at T ~ Z^/^Tq, i.e. at the energy scale of the inner atomic shells. However the effect 
of the deactivation of the inner electrons is small, of order Z~^ in comparison with the total 
number Z of the electrons. On the other hand, at low T, including the most interesting 
region of T ~ Tq, the integral in Eq. (l39|l is determined by the range of x of order one, where 
the model treatment is reasonably justified. 

In order to apply the same model for an estimate of the correlation effect we replace in 
the estimated correlation contribution to the magnetic neutrino scattering in Eq. (l38|) the 
factor |'?/'(0)p by the total density of the electrons that an active electron 'sees' at its location 
in the atom. Then the resulting correction to the integral l\ for an atom can be written in 
terms of the density Uair) of the active electrons and the total density nir) of the electrons 
in the atom: 




(39) 



X = 



2(4:/3nY^^maZ^^^, the energy scale Tq is given by 




(40) 




(41) 



15 



1 



0.95- 

0.9- 
0.85- 

0.8- 
0.75- 

0.7- 
0.65- 

0.6- 

0.55 L 
10 



10" 



10' 



10' 



T (keV) 



FIG. 2: The activation factor / for germanium in the stepping approximation with the actual 
energies of the orbitals (sohd hne) and its interpolation in the Thomas- Fermi model (dashed). 



It should be pointed out that the numerical coefficient in this expression contains a factor 
of 1/8 as compared to Eq.(l38l). This is because of a factor of 1/4 corresponding to the 
statistical weight of the spin-singlet state and an extra 1/2 compensating for the double 
counting of electrons in the pairs. 

One can write the correction described by Eq. (HT]) in terms of a correction to the activation 
factor in the Thomas-Fermi model as 

where the correlation energy scale Ti is given by 

Ti = —ma^ Z ^A.lZeV (43) 

and evaluates to about 131 eV in germanium. The plot of the estimated correlation correction 
in germanium is shown in Fig. 3. One can readily see that this correction is below 1.5% at 



X^/^(x) - ( x{x) - 



T 



X 



3/2 



(42) 
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T ^ 0.3 keV and rapidly decreases at higher energy transfer. Clearly, this estimate refers 
only to the magnetic part of the scattering, while for the weak part we find no correlation 
effect in the considered order due to the cancellation found in Eg. (155]) . We thus conclude 
that in the range of values of T above a few hundred eV the correlation effect can be safely 
neglected for both contributions to the neutrino scattering on germanium. 




T (keV) 



FIG. 3: The correlation correction fc to the activation factor / for germanium in the Thomas-Fermi 
model. 

VI. SUMMARY 

We have considered the scattering of neutrinos on electrons bound in atoms. Our main 
finding is that the differential over the energy transfer cross section given by the free-electron 
formulas ([1]) and (|2]) and the stepping behavior of the activation factor given by Eq. fl2T]) pro- 
vides a very accurate description of the neutrino-impact ionization of a complex atom, such 
as germanium, down to quite low energy transfer. The deviation from this approximation 
due to the onset of the ionization near the threshold is less than 5% (of the height of the 
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step) for the K electrons, if one applies the analytical behavior of this onset that we find for 
the ground state of a hydrogen-like ion. We also find that the free-electron expressions for 
the cross section are not affected by the atomic binding effects in the semiclassical limit and 
for independent electrons. For this reason we expect that the deviation of the actual onset 
from a step function at the threshold for ionization of higher atomic orbitals is even smaller 
than for the ground state, since the motion in the higher states is closer to the semiclassical 
limit. The approximation of independent electrons lacks an account for the two-electron 
correlations arising from the Coulomb interaction between the electrons in the atom. We 
estimate this effect in the large Z limit using the Thomas-Fermi model and argue that the 
effect of the correlations is small in germanium for the values of the energy transfer above 
0.2 0.3 keV. We thus argue that for practical applications, i.e. for the analysis of data 
of the searches for NMM one can safely apply the free-electron formulas and the stepping 
approximation at the energy transfer down to this range. 
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Appendix A: Sum rules 

We consider here the general sum rules for the dynamical structure factor S{T, g^), which 
stem from the analyticity of the density- density Green's function -F(T, q^) at a fixed T and 
complex and also from its asymptotic behavior at large large At a non zero T the 
dynamical structure function, defined by Eq.(|6]), vanishes at = 0, due to the orthogonality 
of the excited states \n) and the initial state |0) in Eq.(|6]) since p(0) reduces to a unit operator. 
For this reason the function F{T, q^) is real at = and thus satisfies in the complex plane 
the condition -F(T, z*) = F*{T, z). At a non zero real q^ the imaginary part of this function 
is not vanishing for both positive and negative g^, so that it has cuts along the real axis 
extending from zero to both infinities {22 1. On the other hand, the asymptotic at large \q^\ 
behavior of the Green's function F{T,q^) is determined by the free-electron formula (|T6|) . 
since at \q^\ — )■ 00 any interaction terms can be neglected. For a scattering on an atom with 
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Z electrons one finds 

This behavior enables one to write a dispersion relation for the Green's function with no 
subtractions: 

By comparing the dispersion relation at — > oo with the asymptotic behavior in Eq. (]Aip 
one readily finds the sum rule for an integral similar to /2, but extended to include also the 
negative c^: 

/oo 
S{T,q^)dq^ = 2mZ , (A3) 
'OO 

where the dynamical structure function at negative is defined by the analytical continu- 
ation and Eq. flTU]) . rather than by Eq.([n]). 

In order to derive from Eq. (lA2p a relation for an integral similar to l\ it is necessary to 
consider the Green's function near the origin, i.e. at — i- 0. In multi-electron systems the 
behavior in this region is generally complicated by the two-electron correlations. For this 
reason we limit the consideration here to the system with just one electron, Z = \. In such 
a system one has p(g) — )■ 1 at g — > 0, so that the Green's function in Eq.(|H]) is contributed 
by only the initial state |0): 

F(T, 0) = i . (A4) 
By comparing this formula with Eq.f lA2p at — )■ one immediately finds the sum rule 

5(T,g2)^ = l. (A5) 
g^ i 

It should be pointed out that unlike the sum rule f lASp this latter relation is generally 
invalidated in multi-electron system by the correlation effects. In fact an indication of such 
a difference in the behavior of the two integrals can be seen in Eq.f l38p . where the discussed 
there correlation effect vanishes for the integral /2, but not for the l\. 

The sum rule (lASP can also be derived from the latter expression in Eq. (fT8|l . Indeed one 
can rewrite the formula as 



oo 



1 (p- g) 



, (A6) 



2m 2m 

and consider the expansion of the last term in powers of (p ■ g). Only the even terms in 
this expansion are non vanishing, since the odd terms give zero due to the parity. One can 
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readily see that in each term in the expansion the pole in is of a higher order than the 
power of in the numerator, so that the imaginary part of each term integrates to zero 
in the integral as in Eq. (lA5p . while the term of the zeroth order in {p ■ q) in Eq. (IA6p gives 
the sum rule flASp . It is again important here that the integration runs over all values of 
i.e. from — oo to +00, since only in this case all the poles of the terms in the expansion are 
within the integration range. Any restriction of the range of integration over may leave 
some poles out so that the vanishing of the contribution of all higher terms in the expansion 
is generally not guaranteed. 



Appendix B: Momentum-transfer integrals for hydrogen-like states 



Consider the situation when the initial electron occupies the discrete nl orbital in a 
Coulomb potential V{r) = —aZ/r. The dynamical structure factor for this hydrogen-like 
system is given by 



S(ni)(T, q^ 



mk 



(27r)3 2/ 



1 ' 

+1 ^ 



mi=-l 



d^k\{vAp{^\Vnlmi)\ 



(Bl) 



where ipnimi is the bound-state wave function, is the outgoing Coulomb wave for the 



ejected electron with momentum fc, and k = \k\ = ■\/2mT — p^, with p„ = aZm/n being 
the electron momentum in the nth Bohr orbit. The closed- form expressions for the bound- 
free transition matrix elements in Eq. ( IBip can be found, for instance, in Ref. 16|. In 
principle, they allow for performing angular integrations in Eq. (IBip analytically. This task, 
however, turns out to be formidable for large values of n. Therefore, below we restrict our 
consideration to the n = 1, 2 states only, which nevertheless is enough for demonstrating the 



validity of the semiclassica' 



Using results of Ref. 



approach developed in Sec. Illli 



17|, we can present the function (IBip when n = 1, 2 as 



S{ni){T, q^) 



2^mp^ 



nl\ 



3[1 - exp(-27rr7)] [{q^ - F + pl)^ + 4p2A;2]2n+i 

2pnk 



X exp 



-21] arctan 



(B2) 



g2 _ J^2 _^p2^ 

where the branch of the arctangent function should be used that lies between and vr. 
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1] = aZm/k is the Sommerfeld parameter, and 

f^^(q^) = 8 [3q'' - {32pl + llA;^)^^ + (82p^ + 72pie + lAe)q^ 
+ {2Qpl - 62pie - 20plk^ - Qk%^ + {pI + k^) 
X {^^pI - ^pte - Iplk' - k'^ + {Apl + e){pl + k'f 

f2p{q') = 2pl [36g« - 48(p^ + k')q' + (152p^ - A%plk' - 8k')q' + {pi + k' 

/1712_^4 ^ 1568_^2l2 , icl4\ „2 , / 44_^2 , /i;^2\ , u2\3 



(B3) 



X 



V 15 



-P2 



^^-pie + lQk' q'+{fpl + Ak'']ipl + e 



(B4) 



(B5) 



Insertion of Eq. ( ]B2p into the integrals (1T5|) and integration over g^, using the change of 
variable 

2pnfc _ ^^^^ 

q^ — k"^ + p2 

and the standard integrals involving the products of the exponential function and the powers 
of sine and cosine functions, yields 



r(l^) 



(T) 



T- 



2mT 



1 - exp( 2^ 
X exp 



-2 



1 — exp 



TT 



: arctan 



1-- + :^)^(B6) 

yi 3yf 



T 



-1 



1 — exp(- 



X exp 



Att 



-4 



1 — exp 



2ti 



Vl/2 - 1 



arctan 



y2-2 



1 

2/2 



80 448 1792 
l%f ^T%| 



(B7) 



r(2s) 



(T) 



2m 



1 — exp(- 



X exp 



47r 



-4 



1 — exp 



27r 



A/y2 



arctan 



\/y2 

2v/y^ 



8 80 

2/2 3^/2 



448 1024 



152/f 



(B8) 



T 



^1 



1 — exp(- 



X exp 



Alt 



-4 



1 — exp 



271 



arctan 



/ 2/2-2 



8 80 
2/2 



704 3328 
+ 



152/i 



451/1 



(B9) 
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i — exp 
X exp 



-4 



V?/2 - 1 



arctan 



v y2 - -L 
( y2-2 



1 

y2 



80 



704 512 

+ 



(BIO) 



where ?/„ = 2mT/p^ = T/lE'nI. The largest deviations of these integrals from the free- 
electron analogs (fT7|) occur at the ionization threshold T = \En\- The corresponding relative 
values in this specific case are 



(Is) 



(FE) 



j{FE) 



1 



7 



0.9572635093, 



j{2s) 
-'1 


= 1 - 


1639 


j(FE) 


15 


r{2p) 
-'l 


= 1 - 


2101 


AFE) 
-'l 


45 



0.9633451168, 



0.9843376226, 



j{2s) 
f2 

■2 

[2 

AFE) 



871 

103 

( 

15 



,-8 



0.9805208034, 
0.9976964900. 



The above results indicate a clear tendency: the larger n and /, the closer I)^'''' and Jg"''' 
are to the free-electron values. The departure from the free-electron behavior does not 
exceed several percent at most. These observations provide a solid base for the semiclassical 
approach of Sec. IIIII 
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